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Abstract 

We discuss and compare two alternative perturbation approaches 
for the calculation of the period of nonlinear systems based on the 
Lindstedt-Poincare technique. As illustrative examples we choose 
one-dimensional anharmonic oscillators and the Van der Pol equa- 
tion. Our results show that each approach is better for just one type 
of model considered here. 
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1 Introduction 



There are several methods for removing the secular terms produced by the 
straightforward application of perturbation theory to nonlinear motion [1]. 
One of them is the Lindstedt-Poincare technique (LPT) [1,2], recently im- 
proved by Amore et al. [3-6] by means of the delta expansion and the principle 
of minimal sensitivity (PMS) [7] giving rise to the LPLDE method. 

There is also another approach that resembles the LPT [8] that we will 
call alternative Lindstedt-Poincare technique (ALPT) from now on. The 
purpose of this letter is the discussion of those alternative approaches and 
comparison of their results for simple nontrivial models. 

In Sec.Hwe briefly review the LPT [1,2] and the LPLDE [3-6]. In Sec. El 
we outline the ALPT [8] . In all those cases we choose the Duffing oscillator [1] 
as an illustrative example. In Sec. 0]we apply the LPLDE and the ALPT to 
other anharmonic oscillators. In Sec. El we apply both methods to the Van 
der Pol equation. Finally, in Sec. El we compare the frequencies of the motion 
for the above mentioned models calculated with both approaches. 



2 The Lindstedt-Poincare Technique 

The LPT consists of the simultaneous expansion of the trajectory and oscil- 
lation frequency in powers of the perturbation parameter. For concreteness 
we choose the Duffing anharmonic oscillator as an illustrative example, [1]: 

fgl + x{ t) = -^(t), (1) 
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with the initial conditions x(0) = 1 and = 0. The parameter fj, is a 
measure of the nonlinearity of the motion. 

Straightforward application of perturbation theory to Eq. (JTJ produces 
an approximate solution in the form of a /i-power series x(t) = Y^^=o x n(t)[i n 
[1,2]. It is well-known that the equations for the coefficients x n (t) have 
resonant terms that after integration give rise to nonperiodic contributions 
[1,2]. Such terms are unbounded in spite of the fact that the solution x(t) 
is known to be periodic for all \i > —1. For this reason, the perturbative 
solution fails at large time scales. 

There are several approaches that provide approximate solutions free from 
secular terms [1]. Here we are interested in the Lindstedt-Poincare technique 
that consists of introducing the actual frequency of motion Q into the equa- 
tion of motion by means of a time dilatation r = Qt [1,2]. Thus the Duffing 
equation (JTJ) becomes 

n 2 X(r)+X(r) = -^(r), (2) 

where X(t) = x(r/0) and the dot represents derivation with respect to r. 
Notice that the period of X(t) is independent of \x. On expanding both X{r) 
and Q 2 in powers of fi, 

oo oo 

n 2 = £ a nf t n , X(t) = £ X n (rK , (3) 

n=0 n=0 

with ocq — 1, one obtains the set of equations 



XJt)+XJt) = -0(n-l) 



j=0 fc=0 



i=i 

n = 0,1,2,... (4) 



where 9 is the Heaviside function. 

If we choose the coefficients a n in order to remove the resonant terms, 
then the solutions X n {r) are periodic, have the general form 

X n (r) = X>n,cos[(2j + I)r], (5) 

3=0 

and satisfy the initial conditions 

X n {0) = <W X n (0) = . (6) 

This approach not only gives us a truly periodic approximate solution but 
also the frequency of the motion in the form of a power series. For example, 
through third order we have: 

flf , = 1 + 3 _ 2 3 , , 

(3) 4 p 12g p 512 p v ) 

Since the motion is unbounded for /i < —1 then the /i-power series is ex- 
pected to be valid only for < I. This is a serious limitation of the LPT 
because the motion is periodic for all /x > — 1 as indicated above. 



2.1 The Lindstedt— Poincare Technique and the Linear 
Delta Expansion 

In order to improve the LPT, Amore et al. [3-6] proposed its combination 
with the linear delta expansion (LDE). The LDE is a variational perturba- 
tion theory like the one proposed some time ago as a renormalization of the 
perturbation series in quantum mechanics [9] which has proved being suitable 
for the treatment of a wide variety of problems [2, 10]. 



4 



In order to apply the LDE to the Duffing oscillator we rewrite Eq. (JIJ) as 
d 2 x(t) 



dt 2 



+ (l + A 2 ) x(t) = 5 [-{ix 3 (t) + X 2 x(t) 



(8) 



where A is a variational parameter and 6 is an order-counting parameter. 
When we set 5 equal to one we recover the original equation (JTJ) that is 
independent of A. Following the LPT outlined above we change the time 
variable and obtain 



Q 2 X(t) + (1 + A 2 )X(r) = 5 -^X 3 (t) + X 2 X{r 



(9) 



Next we expand both X and Q 2 in powers of 5 and proceed exactly as in- 
dicated above for the Lindstedt-Poincare technique, except that in this case 
we choose ao = 1 + A 2 . The resulting perturbation equation of order n reads 



XJt)+XJt) 



(n-1) 



Q'o 



n-1 n—l—j 

E E X 3 {r)X k {r)X n ^ k {r) 

j=0 k=0 



+ Y J ^X n ^{r)-\ 2 X n ^{r) 



n 



0,1,... 



(10) 



and the solutions satisfy the initial conditions As before, the coefficient 
a n is set to remove the resonant term from the equation of order n. At third 
order we find, after setting 6=1, 



3/i 2 ^ 3/i 2 (3^x - 4A 2 



11^ 



J(3) * ' 4 128(1 + A 2 ) ' 512(1 + A 2 ) 2 
We assume that the optimum value of the arbitrary variational parameter A 



is given by the principle of minimal sensitivity (PMS) [7] 



(AT) 



dX 



(12) 



that gives us Xp^s anc ^ an approximate expression of order N. At third 
order we have 



A (3) _ V3^ , , 

A PAfS — 2 ^ ' 



"(3) = On,. > , ■ 



and the approximate frequency results to be 

^ 2 _ 69/i 2 + 192/i + 128 
32(3^ + 4) 

We have already named this approach Lindstedt-Poincare-Linear-Delta- 
Expansion (LPLDE). 

The LPLDE approximants are not /x-power series and prove to be valid 
for all fi > — 1. 



3 Alternative Lindstedt— Poincare Technique 



There is another approach that resembles the LPT discussed above that 
consists of rewriting Eq. (JTJ) as [8] 

d 2 x{t) 

i- I is — 

n=l 



dt 2 



+ In 2 - ]T « n /i n x(t) = -/ix 3 (t) . 



(15) 



Expanding x(t) in a /i-power series and considering f2 independent of this 
parameter, the perturbative equation of order n becomes 
d 2 x n (t) 



dt 2 



+ VL 2 x n {t) 



in- I) 



n— 1 n— 1— j 

I] Z a?i(*)a:*(*)a;n-i-i-fc(*) 

j=0 fc=0 



-f- ^ ] OijX n . 
3=0 



(16) 



The general form of the solutions x n (t) which fulfill the initial conditions 
is 



x n (t) =ZA y cos[(2j + l)fii] 

3=0 



(17) 
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provided that the coefficients a n remove the resonant terms. These coeffi- 
cients turn out to depend upon Q and, as a result, any truncated expansion 
for Q 2 is a self-consistent equation for the frequency. For example, solving 
for Q 2 ^ in the equation of third order 

^- 1 + f-i^ + " 3 ' < 18 > 

we obtain 



n2 _ V30/^ + 96/i + 64 + 2(3/i + 4) 
(3) 16 • 

We have earlier called this approximation ALPT. 



4 Other Anharmonic Potentials 

In order to compare the approaches outlined above we also consider anhar- 
monic oscillators with potentials 

V(x)= 1 -x 2 + ^x 2N , iV = 2, 3,... (20) 

that lead to the equations of motion 



d 2 x(t) 



+ x(t) = -/2x 21 "- L , (21) 



dt 2 

with exactly the same initial conditions considered before for the Duffing 
oscillator, which is a particular case with N = 2. Notice that there is periodic 
motion for all /i > — 1 consistent with the initial conditions. 

Here we apply the LPLDE and the ALPT to the cases N = 3 and N = 4. 
The procedure is straightforward, and we only remark that after removal of 
resonant terms the general form of the LPT and LPLDE solutions is 

X n (r) = 5> n ,cos[(2j + l)r], (22) 

j=0 



whereas for the ALPT we have 

Nn 

x n {t) = ]T Ki cos [{2j + l)Qt) . (23) 

3=0 

5 The Van der Pol equation 

In order to test the performance of the approaches outlined above we also 
choose the well-known Van der Pol (VdP) equation 

^ + x(t) = »[l-x 2 (t)}^. (24) 

because the behaviour of its solutions is completely different from those of 
the anharmonic oscillators discussed above. 

The VdP equation exhibits a limit cycle and leads to oscillations with a 
definite period that depends on the coupling fi. However, unlike the anhar- 
monic oscillators the VdP equation does not correspond to a conservative 
system as the driving term either damps or enhances the oscillation depend- 
ing upon the size of x. This equation has already been treated by means 
of the LPLDE approach [6]; here we apply and compare it with the ALP 
method. 

5.1 The LPLDE Technique 

We first rewrite the VdP equation (|2"2)) as [6] 

Q 2 x{t) + (1 + X 2 )x{t) = 5{fiQ[l - x 2 {t)]x{t) + X 2 x(r)} . (25) 
In this case we expand Q (instead of Q 2 ) in powers of 5 

oo oo 

n=J2°<nS n , x = 5> n (r)(f\ (26) 

n=0 n=0 
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and choose «o = 1 + A 2 . We thus obtain 



xJt) +X n (T) 



(n-1) 
a 



n-1 



3=0 j=l 



n—1 n—l—i n—l—i—j 

"^E E E a i ^(r)x fc (r)i;„_ 1 _ i _ i _ fc (T) 

i=0 j=0 k=0 



.(27) 



The general form of the solutions is 



x n{T~) = E l c nj cos(2j + l)r + d nj sin(2j + l)r] 

3=0 



(28) 



that satisfy the initial conditions 



x n (0) = A n , x n (0) = 



(29) 



for n > 1. The appropriate choice of a n and A n _i enables us to remove all 
resonant terms at order n. 

5.2 The ALPT 

Proceeding as in the case of the anharmonic oscillators we substitute the 
expansion for the frequency Q 2 into the VdP equation and obtain 



d 2 x n {t) 
dt 2 



+ tt 2 x n (t) = 6(n - 1) 



E dx n -\(t) 
UjXn-j{t) + J 

3=1 UL 



jl=0 J2=0 



(30) 



The general form of the solutions in this case is 



x n(t) = E i°nj cos[(2j + + d nj sin[(2j + 

3=0 



(31) 



that satisfy the initial conditions (|29|). 



The appropriate choice of a n and A n ^i enable us to remove secular terms 
at order n. As before, the coefficients a n of the expansion of Q 2 in powers 
of n depend on Q that one obtains as a real solution of the partial sum at a 
given order. 

6 Results and Discussion 

We first consider the results of the LPLDE and ALPT for the Duffing oscilla- 
tor. In the case of the LPLDE we use at all orders the variational parameter 
ApMS calculated at third order. Fig. (JTJ) shows the error over the frequency, 
defined as A = \Q exac t — ^a PP rox\, as a function of // (// > 0) produced by par- 
tial sums through order 20 for both approaches. We observe that the LPLDE 
technique, which yields much more accurate results than the straightforward 
LPT, is less accurate than the ALPT. A similar behavior can be observed 
in Fig. (J2J), where the error over the frequency is plotted for — 1 < // < 0. 
Interestingly, close to \i = — 1, which corresponds to a condition of (unstable) 
equilibrium, the LPLDE method performs better than the ALPT. 

Similar conclusions can be drawn in the case of anharmonic oscillators 
with greater N. Figs. Q and (J3J) show the error over the frequency as a 
function of the order of approximation for quartic (Duffing), sextic and octic 
anharmonic oscillators with \i — 1 and /i = 100, respectively. 

The situation is different in the case of the nonconservative VdP equation 
because the ALPT exhibits a great convergence rate, but it yields a wrong 
frequency for \i >~ 2, which suggests a wrong dependence on this coupling 
parameter. On the other hand, the LPLDE yields reasonable results for any 
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Figure 1: Error over the frequency as a function of \i. 

value of shown in Fig. (jHJ). 

Our results suggest that the ALPT method is more accurate for conserva- 
tive systems but fails for nonconservative ones, while, on the other hand, the 
LPLDE yields reasonable results in both cases. At present we are unable to 
explain the failure of the ALPT for nonconservative systems and we believe 
that this subject deserves further investigation. 
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Figure 3: Error over the frequency as a function of the number of terms 
taking fi = 1 for different anharmonic potentials. 
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Figure 4: Error over the frequency as a function of the number of terms 
taking fi = 100 for different anharmonic potentials. 
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